Coherence and Decoherence in Biological Systems: Principles of Noise Assisted 
Transport and the Origin of Long-lived Coherences 
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The quantum dynamics of transport networks in the presence of noisy environments have recently 
received renewed attention with the discovery of long-lived coherences in different photosynthetic 
complexes. This experimental evidence has raised two fundamental questions: Firstly, what are the 
mechanisms supporting long-lived coherences and secondly, how can we assess the possible functional 
role that the interplay of noise and quantum coherence might play in the seemingly optimal operation 
of biological systems under natural conditions? Here we review recent results, illuminate them at 
the hand of two paradigmatic systems, the Fenna-Matthew-Olson (FMO) complex and the light 
harvesting complex LHII, and present new progress on both questions. In particular we introduce the 
concept of the phonon antennae and discuss the possible microscopic origin or long-lived electronic 
coherences. 



I. BACKGROUND & MOTIVATION 

Quantum mechanics provides the natural lavifs that 
govern the dynamical evolution of atoms and molecules. 
Under well controlled conditions in almost perfectly iso- 
lated systems, quantum coherence and entanglement can 
exist and be manipulated. While there is a clear notion 
of sub-systems when individual units are spatially well 
separated, as is usually the case in systems suitable for 
quantum information processing, this is not necessarily 
so in molecular complexes, where partitioning the total 
systems by means of the choice of appropriate subsys- 
tems is, to some extent, arbitrary. In such situations, 
when the system is static, coherence and entanglement 
can always be made to vanish in a suitably chosen ba- 
sis [I1I2]. Hence, to observe non-trivial quantum effects 
we need to force the system out of its preferred eigen- 
basis, through intervention from the outside, and thus 
probe quantum coherences between eigenstates, as well 
as the quantum coherence properties of their consequent 
dynamics. We can then argue that, while the static sys- 
tem can sustain quantum coherence, it is the dynamics 
generated by the external perturbation which results in 
relevant quantum phenomena and provides the means to 
probe for interesting quantum properties. 

The external perturbation may either be coherent and 
controlled, for example through the action of laser fields, 
or it may be caused by an interaction with an unobserved 
environment. The latter will never be completely avoid- 
able, as perfect isolation of a physical system from its 
environment is impossible. In fact, in bio-molecular com- 
plexes the interaction with uncontrolled environments 
such as vibrations of the surrounding protein matrix, for 
example in the form of thermal fluctuations, is an impor- 
tant driving force of its energetic dynamics and hence its 
functionality. One may therefore expect that the inter- 
play between the internal coherent quantum dynamics of 
the system and the unavoidable presence of noise intro- 
duced by the environment has been optimized by Nature 
during evolution and may be significant for its function. 

Exploring to what extent coherent quantum dynamics 



play a role in composite quantum systems in contact with 
environments, such as bio-molecular complexes and, im- 
portantly, what role, if any, genuine quantum traits may 
play is thus an interesting and timely problem. Explor- 
ing these issues becomes particularly challenging in those 
bio-molecular systems where the strength and timescales 
associated with the environment-induced dynamics dy- 
namics are often comparable to the coherent (Hamilto- 
nian) intra-system dynamics, and cannot therefore be 
treated separately. 

The investigation of these types of questions has gained 
further urgency with the recent discovery of experi- 
mental evidence that excitation energy transport dy- 
namics in the Fenna-Matthew-Olson (FMO) complex 
and other photosynthetic aggregates exhibits surprisingly 
long-lived coherence [3H5]. The FMO complex is an ex- 
ample of a pigment-protein complex, a network through 
which electronic excitations on individual pigments can 
migrate via excitonic couplings. These experiments sug- 
gest the existence of significant quantum coherences be- 
tween multiple pigments through the presence of oscilla- 
tions in the cross-peaks of 2-D spectra which have been 
observed to persist on timescales up to around a picosec- 
ond [S]. This is notable given that the entire excitation 
energy transport in the FMO complex is concluded in 
less than 5ps. 

Theoretical investigations of the role of pure dephasing 
noise in excitation-energy transfer have found that this 
noise mechanism has the ability to enhance both the rate 
and yield of excitation-energy transfer when compared to 
perfectly coherent evolution [HI [7] . An often quoted ar- 
gument states that in the exciton basis typically used 
in previous studies [S], noise-assisted transport is triv- 
ially explained as the result of noise-induced transitions 
between exciton eigenstates, which cause an energetic 
down-hill relaxation towards the reaction centre. While 
this argument can serve to suggest that noise-assisted 
processes might exist in excitation-energy transfer and 
other forms of transport in quantum networks, it falls 
short, as it cannot explain why in some networks noise 
is exclusively detrimental while in others it can support 



transport. Answering such a question requires a more 
detailed study that identifies basic dynamical principles 
of quantum networks in contact with environments. 

Here we would like to provide a fresh look at the 
problem by elucidating the basic mechanisms which give 
rise to noise-assisted transport and identifying where 
noise and coherence respectively play a constructive role 
p, 7, 9 13 . Identifying and understanding these mech- 
anisms at a deep, intuitive level does provide additional 
value even if the individual processes might have been 
known in some previous contexts. Such classifications 
pave the way towards tailored experimental tests of the 
contributions of individual processes [T3]. Furthermore, 
quantum coherence can only be expected to be relevant if 
it is sufficiently long-lived. Present experiments suggest 
lifetimes of coherences between excited states extending 
up to the picosecond range while at the same time the 
coherence between ground and excited states are known 
to last for only about 50-150 fs or so. The origin of 
these long-lived coherences remains to be fully clarified. 
A complete understanding of these fundamental dynam- 
ical principles and the origin of the long-lived coherences 
will allow for the optimization of network architectures in 
order to utilize quantum coherence and noise to enhance 
the performance of artificial nano-structure devices. 



the system (see Ihs of Fig. [2]). Alternatively, one may 
view dephasing noise as arising from the random fluctua- 
tions of energy levels. As a consequence, the fluctuating 
energy levels will occasionally come energetically close 
enough to allow fast excitation energy transfer between 
the sites (see rhs of Fig. I2|. A moderate amount of 
fluctuations serves to enhance the transport while exces- 
sively strong fluctuations of the site energies will make 
the probability smaller for the sites to be energetically 
close. Hence we expect an optimal finite noise strength 
that maximizes transport between two sites. However, 
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II. PRINCIPLES OF NOISE ASSISTED 
TRANSPORT 

In this section we will briefly review the basic dynam- 
ical principles for noise assisted transport following and 
extending the ideas developed in ffj ; lOj fT^Hn| , and fur- 
ther expand them with the new concept of the phonon 
antenna by which coherent dynamics may assist the fine- 
tuning of the network to the spectral density of the fiuc- 
tuating environment. 



A. Basic dynamical principles 



FIG. 1. Left: Local dephasing, for example due to random 
fluctuations of the energy levels, leads to line-broadening and 
hence increased overlap between sites. Right; Viewing these 
fluctuations dynamically, one finds that the energy gap be- 
tween levels varies in time. The resulting nonlinear depen- 
dence of the transfer rate on the energy gap may therefore 
lead to an enhancement of the average transfer rate in the 
presence of dephasing noise. 

the application of excessive noise and the concomitant 
reduction of the effective transition rate between sites 
can also be of advantage as it may effectively block un- 
favorable transfer paths from being followed |12) . 



We have argued that both noise and coherence are im- 
portant for the time evolution of quantum networks in 
the presence of an environment. But, how does the in- 
terplay between these processes occur? 



1. Bridging energy gaps & blocking paths 

Pigment-protein complexes will consist of a number 
of sites whose energies will generally exhibit a certain 
degree of static disorder, i.e. their on-site energies will 
differ from site to site. If this energy difference is larger 
than the intersite hopping matrix element in the rele- 
vant Hamiltonian, then transitions will be strongly sup- 
pressed. Dephasing noise can come to the rescue here as 
it will lead to line broadening thus leading to increased 
overlap between sites without the loss of excitations from 



2. Constructive & destructive interference 

A quantum dynamical system may exhibit a wealth 
of constructive and destructive interference effects. The 
essential nature of this type of effect may be seen in 
the simple network depicted in FigJ2] where the co- 
herent interaction is described by a Hamiltonian, H = 
ELi Ei\'>-){A + ELi •^fc3(|fc)(3| + h.c), where \i) corre- 
sponds to an excitation in site i and we assume J13 = J23. 
An excitation initially prepared in the antisymmetric 
state IV') = (|1) — |2))/'\/2 forms an eigenstate of this 
Hamiltonian which has zero overlap with the site 3 which 
we assume to be coupled dissipatively to a reaction cen- 
ter. Under natural conditions a pigment-protein complex 
is not normally excited in such an antisymmetric state, 
but will tend to receive a single excitation locally, for 
example on site 1. 



The initial state localised on site 1 may of course be 
considered to be an equally weighted coherent superposi- 
tion of the symmetric and the anti-symmetric states, i.e. 
|1) = [(|1) - |2))/y2 + (|1) + |2))/y2]/y2. Thanks to 
constructive interference the symmetric part of the ini- 
tial state will propagate very rapidly into site 3, and from 
there into the reaction center, while the antisymmetric 
part will not evolve at all. Hence the transfer efhciency 
is limited to 50% in this situation. This trapping of pop- 
ulation will be suppressed either via energetic disorder, 
which will release trapped population by inducing coher- 
ent oscillations between symmetric and anti-symmetric 
states, or by environmental dephasing noise which de- 
grades interference effects and therefore destroys coher- 
ent trapping. 



We note that these two mechanisms of noise-assisted 
transport are essentially independent of microscopic 
noises models, and could even be seen as highlighting the 
processes which noise sources must effect in order to drive 
efhcient transport in networks where disorder and/or in- 
terference are important. The precise parameters and 
dynamical evolutions by which this is achieved by an en- 
vironment do of course depend on the noise models and 
simulation techniques employed, but the physical mecha- 
nisms outlined above, which underlie the advantages and 
optimality of noise-assisted transport, are not qualita- 
tively altered by any choice of model or approach. As 
a result, it is reasonable to illustrate the action of these 
mechanisms using fairly simple simulation tools, and in 
section [IlB| we will demonstrate noise assisted transport 
using Lindblad master equations. This approach neglects 
several important factors, such as the effects of finite tem- 
peratures, but for the non-equilibrium conditions we shall 
consider, this approach yields results which both demon- 
strate the dramatic action of noise-assisted transport and 
which are also robust - at the same qualitative level of 
microscopic description as the proposed mechanisms - to 
changes in noise model. 



Realistic systems require much more description, and 
simulations of experimental results certainly require more 
sophisticated models and simulation techniques. Our aim 
in this paper is to slowly build up a picture of the rich 
array of noise-driven processes which can lead to efficient 
energy transport, starting from the most basic ones de- 
scribed above and then incrementally adding additional 
details which lead to new mechanisms. As layers are 
added, so the methods employed will also become more 
involved. A good example of this is given in the next 
section, where we present a new principle of noise as- 
sisted transport which requires a consideration of two 
important properties of system-environment interactions 
which we have so far neglected in this article; the spec- 
tral density of environmental fluctuations and coherent 
vibronic couplings. 




Sink 

FIG. 2. A three-site network in which two sites 1 and 2 are 
each coupled to a third site 3 via an exchange interaction of 
the same strength. Site three is irreversibly connected to a 
sink. In (a) the excitation is delocalized over two sites (red 
and green) with equal probability of being found at either site 
but with a wave function that is antisymmetric with respect 
to the interchange of red and green. This state will not evolve 
due to destructive interference and hence no excitation will 
ever reach the reaction center. In (b) pure dephasing causes 
the loss of phase coherence and the two tunneling amplitudes 
no longer cancel, eventually leading to a complete excitation 
transfer to the sink. 



3. Splitting energy levels - The phonon antenna 

The presence of strong quantum coherent dynamics 
between energy levels can also allow an open system to 
optimize the efficiency of the energy transport. To this 
end, it should be noted that the energy levels of two 
sites that are coupled coherently will split, leading to new 
eigenstates of the global system, one of which is shifted 
upwards and one that is shifted downwards. As a conse- 
quence, one of the dressed energy levels may be moved 
closer to the desired final point of transport such as the 
reaction center. As the other level is moved upwards, a 
further process is required to drive its population towards 
the lower lying energy levels. 

To this end, note that in this new basis dephasing 
noise will now induce transitions between these eigen- 
states (phase noise in the site basis becomes amplitude 
noise in the dressed basis) leading to energy transport to- 
wards the lower lying of the two energy levels. The transi- 
tion rate between these two states will depend on temper- 
ature, the matrix element between these two eigenstates 
and, crucially, the spectral density at the energy differ- 
ence between the two eigenstates. Matching the energy 
level splitting to the maximum of the spectral density of 
the environmental fluctuations can thus optimize energy 
transport. In this sense, we can argue that the two eigen- 
state of the coupled Hamiltonian harvest environmental 
noise to enhance excitation energy transport through the 
formation of a " phonon antenna" . 

As an illustrative example consider again the toy model 
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FIG. 3. In the upper half of the figure, two closely spaced 
energy levels are separated from a third level to which exci- 
tations should be delivered. A coherent interaction between 
the upper two energy levels leads to dressed states |±) with 
an energy splitting which, if matched to the maximum of the 
environment spectral density will optimize transport from the 
upper to the lower level. Hence the split levels act as an an- 
tennae for environmental fluctuations. 




FIG. 4. Population transferred to sink at t = 5ps as a func- 
tion of Ji2 for a Lorentzian spectral density with FWHM of 
50cm~^ (blue dots) and 10cm~^ (purple squares). Numerical 
results obtained using secular Bloch-Redfield equations. For 
all curves J23 — 30cm~^ and the environment temperature is 
77 K. 



consisting in a three site network where site 3 is coupled 
to a reaction center and let us assume that a Lorentzian 
spectral density that is peaked at 200cm^^ (which is typ- 
ical for FMO) and a FWHM of around 50cm^^ (which 
is narrower than estimated for the FMO). Setting the 
energy of site 1 and 2 to ei = £2 = 300c?7i~^ and 
£3 = Ocm~^, it then turns out that the transfer will be 
optimized for a coupling J 12 between sites 1 and 2 of 
around J12 = lOOcm"^ which results in a splitting of 
AE = 200cm^^. Figure |4] shows the population trans- 
ferred to a sink Psink connected to site 3 after 5 ps of 
evolution described by a standard secular Bloch-Redfield 
master equation and the Lorentzian spectral function de- 



scribed above. The transfer to the sink was described by 
a Lindbladian with a transfer rate of lps~^, J13 = 0, 
J23 — 30cm~^ and the bath temperature was taken as 
77K. A clear maximum is seen for J12 — lOOcm"-'^ which 
can be interpreted as a result of the phonon antenna ef- 
fect set out in Fig. [3| Figure |4] also shows results for 
a narrower Lorentzian (FWHM 10cm~^), showing that 
the enhancement at J = lOOcm"^ is much more strongly 
peaked as one would expect. These results suggest that 
exciton couplings might be used to exploit the structure 
of the noise spectrum of the protein surroundings, and 
that the advantages of doing this increase as the structure 
of the environmental spectral function becomes sharper. 

However, when the structure of the spectral function 
contains very sharp resonances, the resonant modes can 
no longer be treated with standard master equation ap- 
proaches. The derivation of most master equations from 
microscopic system-bath Hamiltonians normally relies on 
the fact that the environment contains an extremely large 
number of modes and that the coupling to any individ- 
ual mode is weak. However, a spectral function that is 
very sharply peaked at some characteristic frequency uj 
implies a strong coupling to a discrete mode of the envi- 
ronment, and such an interaction can induce a significant 
coherent modification of the exciton dynamics while the 
mode is driven strongly out of equilibrium by the exciton 
dynamics. These discrete mode effects were explored us- 
ing a recently-developed time-adaptive renormalisation 
group technique that allows strong interactions to dis- 
crete modes to be treated as part of an environment with 
an arbitrary residual spectral density il6 _. This method 
produces numerically exact results with controllable er- 
rors, and it was found that discrete modes may coherently 
modulate population dynamics over ps timescales whilst 
the continuous part of the spectral density drives efficient 
population transfer [THJ [17] . 

Another approach to this problem is to consider the 
strong coupling to a discrete mode as part of the system 
and then to treat the residual smooth parts of the spec- 
tral density as the environment. This approach requires 
a large expansion of the Hilbert space of the effective 
exciton-mode system and makes simulations considerably 
more complex, especially if the residual bath has a non- 
Markovian character. This approach is thus only suitable 
for simple systems with relatively simple residual spectral 
densities, but it is very useful for illustrating the elemen- 
tary concepts of strong exciton-mode interactions. We 
now consider such a model using the Hamiltonian given 
by Eq. ([2) to describe three sites, as described above, in 
which each site is coupled coherently to a single harmonic 
oscillator of frequency lu — 200cm^^ and coupling con- 
stant gi — 30cm^. To clearly demonstrate the discrete 
phonon antenna effect, we neglect background dephasing 
here, but this will be treated in detail in a forthcoming 
work |18) . The simulations are performed using exact 
diagonalisation and the simulation results were found to 
converge with 4 bosonic levels per mode. 

Figure [5] shows the population in the sink as a func- 
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FIG. 5. Paink as a function of time for three-site system, 
where each site is coupled to a mode of frequency 200cm~^ 
with coupling strength g = 30cni~^. A strong enhance- 
ment of transport is seen for J12 = lOOcm"^ (squares). The 
other curves are J12 = 150cm~^ (diamonds), J12 ~ 200cm~'^ 
(triangles), Ji 2 — SOOcm"^ (inverted squares) and J12 = 
50cm^^ (dots). For all curves J23 — 30cm~^ and the initial 
state of the vibrational mode is a thermal state at 77 K. 
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FIG. 6. Psink as a function of time for a three-site system 
without any coupling to a mode. Weak transport is now 
seen for J12 = lOOcm"^ (squares), the other curves are 
J12 = 150cm~^ (diamonds), J12 = 200cm~^ (triangles), J12 = 
300cm~^ (inverted squares) and J12 = 50cm~^ (dots). For all 
curves J23 = 30cm~^ and the initial state of the vibrational 
mode is a thermal state at 77 K. 



tion of time for a range of J12 values. It can be seen 
that the behaviour of Psink is strongly non-nionotonic 
with a large enhancement of the population transport 
when J12 = lOOcm^^. For these values the hybrid states 
|±) are resonant with the vibration sidebands of each 
site (see Fig. [7| and the mode-exciton interaction medi- 
ates strong transitions between sites by allowing efficient 
coherent tunneling between the degenerate vibrationally 
excited states of the electronic excitations. For other 
values of J12 the |±) states are detuned from the mode 
sidebands and vibration-mediated transport is dramat- 
ically reduced. Note also that for J12 — 300cm~^ the 
lower state |— ) is resonant with site 3. This leads to a 
rapid transport of 50% of the population, but the 50% 
initially held in the |-|-) state is at such a high energy 
that it cannot transfer population to the sink via site 3 



Il> 



\2> 



2J 



|-> 



* * 



sink 



1|3= 



FIG. 7. Energy level scheme for three-site and sink system 
(blue lines) showing vibrational sideband levels (dashed red 
lines). The coupling J12 splits the degenerate transitions, 
forming new super position states |±), and can bring these 
electronic levels into resonance with the sidband structure. 
Under these conditions fast, coherent inter-site transport 
(black arrows) becomes possible through the new exciton- 
mode energy level landscape. 



in the absence of background noise. 

Figure |6] shows results for the same parameters in the 
absence of the mode. One now observes a monotonic be- 
haviour of the transport rates as J12 is increased, and 
by comparison with Fig. [5] the dramatic resonant en- 
hancement of transport observed at J12 = lOOcm"^ is 
highlighted. These results suggest that tuning of exci- 
tonic couplings to "harvest" the available noise spectrum 
is most effective when the quasi-discrete features are ex- 
ploited. It is in this sense - just as excitonic interactions 
allow chromophores to sample larger parts of the opti- 
cal spectrum - that we refer to such coherent sampling 
of discrete features or maxima of the spectral function 
as the phonon antennae effect. Not only do these allow 
for tremendous gains in transport efficiency, their strong 
quantum mechanical interaction with the exciton appears 
to strongly modify the behaviour and lifetime of coher- 
ence dynamics (not shown). 



B. The principles at work: FMO and LHII 

Here we will briefly analyse how these principles are 
put to work to explain aspects of noise assisted transport 
in two different pigment-protein complexes. 



1. The FMO complex 

To illustrate how the general principles outlined in the 
previous section may feature in a realistic scenario, we 
will consider the example of excitation transfer across a 
monomer of the FMO complex, which can be modeled 
as a 7 sites network and for which detailed information 



concerning the system Hamiltonian is available [TS]. We 
ignore here the recent addition of an eigth site [20] as it 
does not affect the general arguments and was unlikely 
to have been present in recent experiments reported in 

[3|S]. 

Given the strong coupling of sites 1 and 2, these levels 
are shifted and mixed and the dynamics is conveniently 
described using an hybrid basis [1.2^ for these two sites 
that we denote by |±). In this basis, the Hamiltonian 
has the local site energies and coupling structure shown 
in FigurelSk, where site 3 is connected to a sink node from 
which excitation is transferred to the reaction center and 
all the remaining sites have been packed in a block (la- 
beled as "additional sites" in Figure 2) that is uncoupled 
from the level |+). 

When an excitation is injected in site 1, the coherent 
evolution leads to transfer efficiency Pgink below 60%, far 
from the ideal transfer, represented in Fig. |8]by the red 
horizontal line. The coherent interaction between sites 
1 and 2 leads to level splitting and moves one of these 
energetically closer to site 3 while the other is farther 
removed. The introduction of dephasing noise with a 
given spectral density can dramatically add to this pic- 
ture and lead to perfectly efficient transport to the sink, 
as illustrated in the part (b) of Fig. Is] Environmental 
fluctuations lead to transitions between states \+) and 
|— ) and the resulting transition rate will depend on the 
spectral density at the level splitting of the dressed states 
|-|-) and |— ). Indeed, inspection of the relevant Hamil- 
tonian reveals that the coupling between sites 1 and 2 
is of the order of 93cm~^ and leads to a splitting be- 
tween the dressed states |-|-) and |— ) of very close to 
200cr7i~^ which matches the maximum of the spectral 
density of the environmental fluctuations used in 19 . It 
is hence tempting to speculate that nature has optimized 
the intersite coupling and/or spectral density to provide 
a perfect match between the maximum of the spectral 
density and the level splitting of the dressed states |-|-) 
and |— ) thus exploiting the idea that quantum coherent 
interaction support an "antenna for phonons". 

Furthermore, coherent oscillations between level |— ) 
and the rest of the complex are now largely suppressed 
(due to noise and destructive interference), while an inco- 
herent transfer path between previously decoupled levels 
|±) is now active, leading to a fast transfer of popula- 
tion to the sink via |— ). Within this model, it is possi- 
ble to optimize the local dephasing rates so as to repro- 
duce the observed transfer rates in the correct time scale. 
More detailed discussions can be found in the literature 

(aEniniiia. 

Recent molecular dynamics simulations have also 
shown that the spectral density of pigment-protein com- 
plexes such as the FMO complex contain a large num- 
ber of sharp features which could be exploited via the 
mechanism described above [2TH23] . Moreover, spectral 
hole burning and fluorescence line narrowing experiments 
have also shown that a number of modes with frequencies 
comparable to the energy differences between excitons 



couple strongly to the electronic degrees of freedom in the 
FMO complex [24VI26J . Given that vibration-mediated 
tunneling is also implicated in a range of other impor- 
tant biological transport processes, such as olfaction |27| . 
it does not seem unreasonable that EFT dynamics in 
pigment-protein complexes could take advantage of the 
rich spectral structure of their local environments by ap- 
propriately tuning the inter-site interactions to sample 
the maxima of the spectral function. 

Indeed, it has also been shown in [12] that a discrete 
mode coupled to each site in the FMO complex could 
also induce cfhcicnt transport, and in some cases could 
even outperform the optimised, purely incoherent noise 
assisted transport which arises from a Lindblad treat- 
ment of the environment. Given the observed presence 
of a wide range of discrete features in the phonon spec- 
trum of FMO, it is likely that both incoherent and dis- 
crete phonon antenna effects could play a role in the effi- 
cient energy transport seen in experiments, with the lat- 
ter mechanism possibly contributing to enhanced, longer- 
lasting coherent dynamics, as observed in [12]. A full 
analysis of the dynamics of FMO in the presence of both 
discrete modes and continuous spectral densities with 
maxima close to exciton energies differences will be ap- 
pear in [28) . 



2. The LHII complex 

Here we consider the LHII complex from the photosyn- 
thetic unit of purple bacteria [29]. The photosynthetic 
unit is composed of a membrane that contains only two 
types of antennae complexes, LHI and LHII. The former 
are believed to form ring-like structures that enclose the 
reaction center while the latter transfer energy to the 
reaction center via the LHI antennae. The LHII anten- 
nae consists of two subunits, the B80G and the B850 ring. 
The B80G ring, absorbing around 800nm, consists of nine 
pigments, which have their molecular plane arranged per- 
pendicular to the symmetry axis. The B850 ring, absorb- 
ing around 850nm, consists of nine repeating pairs of a-/3 
pigments, with opposite transition dipole moments and 
their molecular planes parallel to the symmetry axis. 

Here we make the simplifying assumption that the in- 
teraction between pigments is well approximated by a 
dipole-dipole interaction of the form 



^'' 4^e|r;jf 



fii^] ~ S{i^irij){fiirij) 



(1) 



where fij — fi — fj is the vector joining pigment i 
with pigment j and fij is the corresponding unit vec- 
tor. We chose e = 1.3 • SMIQ-^'^C^/Nm^. The dipole 
moments are assumed to be of equal modulus |/i| = 
2.046 X lO*-"^^^ Cm. The dipole moments are nearly per- 
pendicular to the symmetry axis and to the radial direc- 
tion (see [30] for details). 

With this information it is possible to determine the 
interaction strengths between any pair of pigments. In 
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FIG. 8. Noise-assisted energy transfer across FMO can be 
qualitatively understood by introducing a hybrid basis of local 
sites. While a purely unitary evolution yields to inefficient 
transport, as discussed in the text, the presence of dephasing 
noise eliminates inefficient transport paths while opening up 
new channels for excitation transfer. While the estimation for 
the transfer efficiency under coherent evolution is well below 
fOO%, a simple noise model brings this number close to perfect 
efficiency within the observed transfer time. 



the B800 ring neighboring pigments experience an inter- 
action strength of around —26.7cm^^ while in the B850 
ring neighboring pigments within a a-/3 pair experience 
an interaction of about 190.8c?7i~^ and between neigh- 
boring pairs of 239. 4cm^^ while the couphng from B800 
to B850 between the closest pigments is of the order of 
12c77i~^. Employing these parameters it is now straight- 
forward to observe that the transfer from the B800 to the 
B850 ring benefits from the presence of dephasing noise. 
Fig [9] shows the non-nionotonic dependence of the trans- 
fer rate between the two rings versus the noise strength 
under the assumption that all pigments are subject to 
noise of identical strength. Here the noise assists in bridg- 
ing the energy gap between the B800 and the B850 ring. 
Strong noise does however tend to suppress delocaliza- 
tion of excitons in each ring and thus suppress coherent 
superpositions that may lead to enhanced transport from 
B800 to B850. 




FIG. 9. The B800-B850 transfer rate versus the noise strength 
exhibiting a non-monotonic behaviour in the noise strength. 



III. REFLECTIONS ON THE MICROSCOPIC 
ORIGIN OF LONG-LIVED COHERENCE 

Having given an overview of the key concepts underpin- 
ning noise- assisted transport, we proceed with our gen- 
eral discussion of the role of noise in PPC dynamics by 
going a little more into the details of the microscopic 
couplings which generate the essential processes of pure 
dephasing (destruction of interferences) and exciton re- 
laxation. Without leaping into major simulations, we 
aim to point out a number of effects which might be ob- 
servable in experiments, and which could potentially help 
us to pinpoint the physics behind a number of currently 
ill-understood observations that have been made in ex- 
periments and ab initio simulations. 

Let us consider the following standard Hamiltonian of 
chromophores which are linearly coupled to protein fluc- 
tuations in their excited states. Defining bosonic annihi- 
lation and creation operators aki, a{,j for the independent, 
uncorrelated fluctuations of frequency a;^ of the protein 
environments acting on site i, and denoting an optical 
excitation on site i by the state |z), the Hamiltonian for 
a single excitation can be written as. 



H = Y,{{^r+Xim{i\+Hi 



T.J^A^){jl (2) 

where X, = Y.k9ik(.aki + a^i) and H^ = Sfe'^fcoL^fci' 
where ujk are the frequencies of the protein fluctuations 
and gik their coupling strength to state |i) {h = 1 
throughout). Assuming identical environments for each 
site, the dissipative effect of the environment on the ex- 
citons is fully-characterised by the spectral density of the 
environment G{ui) = X]fe5'?fe'^('^ ~ ^k): which we dis- 
cussed extensively in the previous section. It is nor- 
mally taken that the electronic couplings J^- between 



sites is stronger than the pigment-protein interaction and 
that the appropriate basis for analysing experimental 
results corresponds to that of the excitonic eigenstates 
\eu)=Y.rC:\i). 

Rewriting the Hamiltonian in the exciton basis, we ob- 
tain, 

H = ^£^n|e„)(e„| +^Q„|e„)(e„| 

n n 



E^/ 



(3) 



The new couplings to the bath are given by Qn = 

T,k,^9ik\Cl'\^iak^ + al^) andRnm = T,k,i9zkC^C^iakz + 

a|,j). As was mentioned in the introduction, in this ba- 
sis the exciton states, which is the appropriate basis for 
weak system environment coupling, are stationary under 
the action of the electronic Hamiltonian and dynamics 
arise solely through the action of the bath couplings on 
the system. From Eq. pi one can see that these interac- 
tions come in two distinct flavours, namely inter-exciton 
(transverse) transfer terms (ex Rnm) and diagonal (lon- 
gitudinal) fluctuations of the exciton energies (ex Q„). 
During the evolution of an initially-prepared superposi- 
tion of exciton states, the total dephasing time of the su- 
perposition roughly contains two contributions, a relax- 
ation part coming from transverse processes and a pure 
dephasing part coming from the longitudinal processes. 
Before going on to discuss the case of uncorrelated fluc- 
tuations, we briefly present the parameters 



Inm = 2^0, Gj Cj Cj e ''' ", 

- 2(Cf )^(C")^)e-'''^/''% 



(4) 



(5) 



where Vc is the correlation length of the fluctuations at 
each site and Vij are the inter-pigment distances |19j . 
The off-diagonal quantities jnm are proportional to the 
strength of the transition rates between excitons n and 
m in Bloch-Redfield theory, and pure-dephasing pro- 
cesses between excitons n and m in modified Redfield 
theory are proportional to a„m |19j . Although we will 
not discuss spatial correlations explicitly, we note that 
when Tc is much larger than the size of the complex 
Jnm ~ |(en|em)P, a,nm ~ and the total (transverse and 
longitudinal) coupling to the bath vanishes in the one- 
excitation manifold. Spatial correlations thus reduce de- 
phasing and transport rates by roughly similar amounts. 
Returning to uncorrelated environments, it is well 
known that delocalisation of the excitons over extended 
regions of roughly N sites leads to a 1/iV suppression 
of the pure-dephasing of ground state-excited state co- 
herences via the phenomena of motional narrowing |31j . 



However, due to disorder the delocalisation length in typ- 
ical PPCs is often restricted to only a few sites, and in- 
deed several PPCs only contain a few sites in total. Nev- 
ertheless, even if the excitons are only delocalised over 
two sites, a significant suppression of zero-quantum co- 
herence dephasing can arise if the electronic wavefunc- 
tions are nearly equally distributed over the two sites, 
i.e. the eigenstates are symmetric and anti-symmetric 
combination of the two local excitation states. To see 
this clearly, imagine we have a dimer system in which 
ei — £2 ^ J. The two excitonic states are then |e±) = 
(|1) ± |2))/-\/2. In the exciton basis the exciton-phonon 
Hamiltonian becomes. 



+ (i?+_|e+)(e_|+h.c) + ^iJ^ 



|e+)(e+K|e_)(e_| 



=i.d 



(6) 



As indicated by the underbrace, the longitudinal coupling 
to the excitons vanishes (in the one-excitation sector) 
when the eigenstates are symmetric/anti-symmmetric 
combinations. The bath-induced dynamics now pro- 
ceeds purely through transverse terms and dephasing 
only arises through relaxation. We shall refer to this 
case as relaxation-limited dephasing. In general, if the 
excitons are delocalised over N sites and excitonic cou- 
plings are very strong, then the eigenstates are of the 

form |e„) = X]a:=i ^~*"'^^^l*) ^^'^ "^i^l ^^^^ ^^^ tie sub- 
ject to any longitudinal bath coupling. 

Our discussion in this section has so far made no as- 
sumptions about the spectral function of the environ- 
ment, but in order to illustrate a few ideas we will now 
consider the case of a sufficiently weak exciton-phonon 
coupling where the influence of the bath can be treated 
in a Markovian (Redfield) approximation. In this limit, 
one can rigourously define the relaxation time Ti, the 
total dephasing time of coherences T2 and the pure- 
dephasing time T2*. The three times are related by 
I/T2 = l/2Ti + l/T^. The relaxation time Tf^ de- 
scribes the rate of population transfer induced by the 
transverse couplings in Eq. (l3| and T2* arises from lon- 
gitudinal terms. In the limit of relaxation-limited de- 
phasing, where T| diverges, T2 = 2Ti and the dephasing 
rate is actually slower than the population transfer rate. 
However, in most experiments performed in liquids and 
the solid state, T2 <C Ti due to the dominant contribution 
of pure dephasing terms. This is often also true in real- 
isations of qubits in quantum information science, and 
removal of these processes is currently a topic of major 
research. 

In current 2D spectroscopy experiments it is often pos- 
sible to meaure the population lifetimes (giving an esti- 
mate of Ti) and the total dephasing time (T2) of co- 
herences seperately. By comparing these numbers, the 
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extent of pure dephasing can be inferred and the ex- 
tent of the suppression of pure-dephasing terms by de- 
locaHsation might also be extracted. This might possi- 
bly allow strongly-coupled excitonic states to be identi- 
fied and could also thereby provide a consistency-check 
of Hamiltonians and dynamical models which have been 
proposed for these systems. From this point of view, it 
must be noted that non-Markovian dynamics are thought 
to play an important role in PPCs such as FMO [5^ . 
and in these cases there is no simple relation between 
the typical timescales associated with pure dephasing, re- 
laxation and the total dephasing timescale(s). However, 
non-perturbative simulation techniques can address this, 
and this will be dealt with in a forthcoming work [28] . 

One example of this effect could be found by looking 
at excitons associated with sites 1 and 2 of FMO. In 
most published Hamiltonians the difference between lo- 
cal excitation energies of sites 1 and 2 are smaller than 
the large excitonic coupling between then, and two ex- 
citon states contain strong contributions from symmet- 
ric and anti-symmetric combinations of sites 1 and 2. 
Within the model and modified Bloch-Redfield approx- 
imation presented by Adophs and Renger [TS], the co- 
herence pnm between excitons n and m decoheres un- 
der pure dephasing processes as pnm oc 6"°'""*^^^*^"^^°^^ 
( we have neglected the relaxation-induced contribu- 
tion to the decoherence) . The coherence between ex- 
citon n and the ground state pg„ similarly decoheres 
as pgn oc e"'''""^^^*)"^^''^^. In both cases, the spectral 
density-dependent functions g{t) are the same [T^. The 
relative strengths of the pure dephasing dynamics of the 
inter-exciton and ground-exciton can be therefore be es- 
timated by looking at the ratios anm/lnn or anm/lmm, 
which are independent of the spectral function. Lower 
and upper estimates of these ratios were obtained by 
computing a„m,7„„ and 7mm for the excitons |e„), |em) 
of FMO Hamiltonians with the largest contributions from 
sites 1 and 2. Assuming uncorrelated environments, we 

found that anmhnn ~ anm/lmm « 10"^ - 10~^, with 

the lower bound being obtained from the Hamiltonian of 
Adolphs and Renger |19j and the upper bound from the 
Hamiltonian presented in [33] ■ The large suppression of 
the inter-exciton pure dephasing contribution at the up- 
per bound, would most likely lead to a relaxation-limited 
dephasing. 

After completion of this work, an explicit example 
of relaxation-limited dephasing has also been given in 
[34] , They have computed the dynamics for site popu- 
lation dynamics of sites 1 and 2 of the FMO complex 
using a weak coupling non-interaction blip approxima- 
tion (NIBA) developed in the context of the spin-boson 
problem '35] . Pachon and Brumer [34' point out for the 
parameters of these sites, the small site energy difference 
and large excitonic coupling leads to a suppression of 
pure-dephasing contributions to the total dephasing rate, 
and at 77 K, the dephasing rate they predict is indeed 
close to the relaxation-limited case T2 ~ 2Ti. However, 
as we stress above, this relaxation-limiting suppression of 



pure-dephasing - which does not rely on a non-Markovian 
treatment of the system-bath interactions, is only rele- 
vant under conditions where site energies differences are 
much smaller than inter-site couplings. In the FMO com- 
plex experiments, the long-lasting excitonic coherences 
are observed between the two lowest energy excitons, 
which correspond to sites 3 and 4, rather than sites 1 
and 2. These sites have energy differences and couplings 
of roughly equal magnitude [H] , or much greater site en- 
ergy differences [33]. While relaxation-limited dephasing 
can play a role for certain pathways in the FMO dynam- 
ics, it is unlikely to provide a general explanation for 
the most significant long-lasting coherences observed in 
FMO, as is claimed in [34]. It should also be pointed out 
that the long-lasting inter-exciton (not to be confused 
with inter-site) coherences at 77 K in FMO persist for 
approximately 1.8 ps, much longer than those predicted 
in iij. 

It is important to also note that the long-lasting co- 
herences found in |34j are not due to any non-markovian 
effects; the use of the NIBA approach (which may not be 
reliable in the relevant parameter regime, see footnote 
[55] gives non-perturbative renormalisation of the pa- 
rameters which determine the relaxation and dephasing 
rates, but the dynamics themselves are simple, strictly 
exponentially-damped oscillations and the coherence life- 
times could be predicted from polaron perturbation the- 
ory. This, and the agreement of this theory with numer- 
ical simulations, could also suggest that for the simple 
Ohmic spectral density and weak-coupling parameters 
used in the literature, advanced simulation techniques 
are not required for complexes such as FMO, though a 
comprehensive understanding of EET dynamics is still 
lacking as it is unlikely that the simple Ohmic spectral 
density describes all aspects of the system-bath physics. 

One of our reasons for discussing dclocalisation- 
induced suppression of pure-dephasing is that the in- 
terexciton coherences in FMO appear to dephase much 
more slowly (~ 1 ps) than the typical (^ 70 — 150 fs) de- 
phasing times of the ground-exciton coherences between 
excited and ground states [5] . If the dephasing of ground- 
exciton coherences comes from the same set of fluctua- 
tions, one would expect roughly similar pure-dephasing 
times for the zero-quantum coherences as well. One po- 
tential solution is that the fluctuations are spatially corre- 
lated, and has been used to analyse the results of a num- 
ber of experiments on PPCs and conjugated polymers 
|37[l38j. However, a number of recent quantum chemistry 
and molecular dynamics simulations of the FMO protein 
have, somewhat inconveniently, found no evidence for sig- 
niflcant spatial correlations between different pigments 
[2TH23]. 

An alternative possibility that may be appropriate 
for excitons which are strongly delocalised over a few 
sites in PPCs is suggested by relaxation-limited dephas- 
ing. We have shown that strong excitonic coupling can 
suppress pure dephasing amongst excitons in the one- 
excitation sector, but when the electronic ground state 
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\g) is also included in the Hamiltonian, it can be seen 
that the excitons are no longer protected from the lon- 
gitudinal terms if they are in a state that spans exci- 
tation number sectors e.g. a ground-exciton coherence 
|e„)((7|. This means that a ground-exciton coherence 
will feel a motionally-narrowed, but still potentially large 
pure dephasing, whereas an inter-exciton coherence be- 
tween |e±)-like states will not be subject to this pure 
dephasing at all. 

It is important to notice that, unlike spatial correla- 
tions which lead to a dynamical decoupling of the entire 
exciton-phonon interaction, the relaxation-limited case 
corresponds to the entire weight of the exciton-phonon 
coupling being transferred to the inter-exciton couplings. 
In this sense, and under appropriate conditions, coher- 
ences are protected without effecting transport efficiency. 
The inter-exciton coherence dephasing time is now essen- 
tially determined by the lifetimes of excitons, and has a 
much stronger, and thus tunable, dependence on the exci- 
ton energy differences compared to pure dephasing which 
is frequency-independent. For example, the Redfield re- 
laxation rate r_| between excitons -|- and — is given by 

m, 



difference depending sensitively on the spectral density 
used. 

The mechanism described above is another example of 
how excitonic interactions can be used to shift the domi- 
nant effects of the environment between spectral features 
- here removing effects of aggressive low frequency de- 
phasing fluctuations - and conceptually allowing lifetimes 
and coherence times on the same order of magnitude as 
each other and the transport time (which is controlled by 
the lifetimes). This could be an important part, though 
certainly not the only mechanism, of the physics under- 
lying long-lasting cohcrcrnces. However, the question of 
why these systems might choose to do these things - espe- 
cially the preservation of coherences - remains as elusive 
as ever. 

We finally note that this idea of electronic protection of 
inter-exciton coherences could be tested in dimer systems 
with strong excitonic interactions such as those found 
in reaction centers, and a detailed theory of the ideas 
presented in the last two section of this article will appear 
shortly [28]. 



IV. CONCLUSIONS 



r+_ = ^[(1 + n{E+^)G{E+^) + n{E^+)G{E^+)], (7) 

where Enm — En — Em, En is the energy of the nth ex- 
citon and n{uj) is the thermal bosonic occupation of a 
mode with frequency uj. The total dephasing time of an 
inter-exciton coherence can be significantly longer than 
the ground-exciton dephasing time in the relaxation- 
limited case when the exciton lifetimes determined from 
Eq. (l7| are smaller than the ground-exciton dephasing 
rate. For our dimer example, when SE± « 2J is much 
larger than the maximum of G{lu), the lifetimes of the 
states will be long because r„„i ex G{Enm)- The es- 
sential point here is that relaxation and pure-dephasing 
are determined by very different regions and properties 
of the spectral function, with pure dephasing rates at 
long times being very sensitive to the functional form 
of lim^_i.o G(aj) coth(/3a;/2)/w^ [39. Therefore for tran- 
sitions between states separated by large energy differ- 
ences, relaxation-limited dephasing may significantly in- 
crease the effective inter-exciton dephasing time com- 
pared to the ground-exciton dephasing rates, this the 



In this work we have summarized recent work concern- 
ing several aspects of the dynamics of quantum networks 
in contact with environments such as pigment-protein 
complexes. We have discussed the basic principles under- 
lying the dynamics of such systems and formulated the 
concept of efficient transport mediated by phonon anten- 
naes induced by coherent dynamics. We also discussed 
the possible relevance of motional narrowing on the un- 
expectedly long lifetime of zero coherences in pigment- 
protein complexes such as FMO. A more detailed account 
of these two aspects will be pubhshed elsewhere. 
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